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STATUS OF FLUTTER OF FLAT AND CURVED PANELS 
By Robert W. Leonard and John M. Hedgepeth 


SUMMARY 


Representative results are presented to show the current status of 
the panel flutter problem. The discussion includes flat panels with 
and without midplane stresses, buckled panels, and both unstiffened and 
stiffened infinitely long circular cylinders. 


INTRODUCTION 


The flutter of the skin panels of supersonic airplanes and missiles 
has aroused considerable interest in recent years. Numerous investiga- 
tions (refs. 1 to 28 ) have determined that panel flutter may govern the 
design of at least very thin panels such as those used in fairings and in 
radiation shielding. Representative results of some of these investiga- 
tions, which tend to illustrate the current status of the problem, are 
presented in this report. The discussion is divided roughly into three 
parts: First, the flutter of flat panels; second, the flutter of panels 

stressed beyond the buckling load; and third, the flutter of thin-walled 
circular cylinders . 


SYMBOLS 


Ap - 




( 2 ) 


a panel dimension in chordwise direction 



af V_ 

T (2) 

-Lrro 


b 
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panel dimension in spanwise direction 


panel stiffness, 
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Young's modulus 


^J (I) = fo ^ (l) 5?fc (s) ] d5 


1 > J > m >V 
M 

n 


R = St 


integers 
Mach number 

integer, number of circumferential waves around the cylinder 
chordwise midplane load per unit spanwise length 
spanwise midplane load per unit chordwise length 
dynamic pressure 


r 

t 

w 

X 

y 

A 


1(e) 


cylinder radius 
thickness 
panel deflection 
chordwise coordinate 
spanwise coordinate 
dynamic-pressure parameter. 


2qa' 


7m 2 - 1 D 

Poisson's ratio (taken equal to 1/3 throughout) 

mth mode of uniform clamped-clamped beam of unit length 

dimensionless coordinate 


a? 


+ 2 


dx^y 2 c)y 4 
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FLAT PANELS 


Up to the present time, the largest research effort has teen devoted 
to theoretical study of flat panels, Typical analytical results for flat 
panels are shown in figures 1 and 2. 

The results in figure 1 apply to panels of infinite aspect ratio, 
that is, panels which extend to infinity in the spanwise direction and, 
hence, behave structurally like a beam. Figure 2 contains results for 
a panel with finite aspect ratio. Both plots give the panel thickness 
ratio t/a necessary to prevent flutter, at various Mach numbers, of 
steel panels at sea level. The results are restricted to panels which 
are unstressed; that is, they have no midplane tension or compression. 

In both figure 1 and figure 2, the dashed curves at low supersonic 
Mach numbers are two-mode results based on the use of unsteady linearized 
air forces. (See refs. 8, 10, and 20.) The solid curves at high Mach 
numbers are two-mode results obtained by using much simpler "static" air 
forces in which no account is taken of unsteady effects. (See refs. 22, 
26, and 28.) Experience has shown that, for the high Mach numbers, the 
static approximation yields results in close agreement with those yielded 
by more refined theories. The smooth merging of the dashed and solid 
lines for each configuration illustrates this close agreement. 

In figure 1, curves are shown for three different configurations: 
an array of panels continuous over equally spaced supports, a single 
panel with pinned ends, and a single panel with clamped ends. The general 
effect of panel boundary conditions may be seen from the relative posi- 
tions of the curves at high Mach numbers. The most critical configura- 
tion, requiring thicker panels to prevent flutter, is the array of panels 
on equally spaced supports. 

The effect of aspect ratio on required panel thickness is shown in 
figure 2. In this case, both curves are for a single panel with pinned 
edges. The upper curve applies to an infinite span panel and is the 
same as the middle curve of figure 1. The lower curve applies to a 
square panel. As would be expected, reducing the aspect ratio from 
infinity to one reduces the required thickness. 

It should be pointed out that the moderate reduction in thickness 
ratio at high Mach numbers is due primarily to the structural effects of 
aspect ratio. (See ref. 26.) The aerodynamic effect is less than 2 per- 
cent throughout the range of the solid curves. Thus, the use of aero- 
dynamic strip theory would give accurate boundaries for the square plate 
for these high Mach numbers. 
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The situation is apparently different for the low Mach numbers. 

In addition to the moderate structural effect, there is now a large 
aerodynamic effect of aspect ratio and a resulting large reduction in 
thickness in changing from the infinite-aspect-ratio panel to the square 
panel. It is interesting to note that, although the low supersonic Mach 
number range appears to be critical for infinite -aspect-ratio panels, 
this may not be the case for the simply supported square panel for which 
the results show a steady increase in thickness with increase of Mach 
number . 


It should be pointed out that the dashed curves at the low Mach 
numbers in figures 1 and 2 are based on a rather limited number of cal- 
culations and this region has not yet been carefully explored. It is 
known, however, that structural damping is fairly effective in reducing 
the large thicknesses required for infinite -aspect-ratio panels. (See 
ref. 10.) 

It is also worthwhile to note that a small amount of experimental 
data has been reported at Mach numbers of 1.3 and I .56 for clamped panels 
which behave like the infinite-aspect-ratio panel. (See refs. 6 and 13.) 
The data are not shown in figure 1 because they apply to a different 
altitude and material. When compared with calculations made on the same 
basis, experiment and theory agree at Mach number I. 56 . At Mach num- 
ber 1 . 3 , the data show an increase of thickness but not the large increase 
predicted by the theory. The difference may be due in part to structural 
damping. 

So far, only panels with no midplane stress /have been considered. 
However, the effect of midplane tension or compression stresses on the 
required thickness of flat panels has been investigated and is found to 
be important (ref. 26) . This effect is illustrated by the results shown 
in figure 3* 


Figure 3 is a plot of a modified- thickness-ratio parameter 

containing both panel and air properties, against a 


chordwise compression parameter 


(Buckling P x )p = q 

y 


The coefficient of 


the thickness ratio contains the Mach number M, Young's modulus E, 
and the dynamic pressure q and allows the application of the plotted 
results to all Mach numbers above about 1.5* The compression parameter, 
which specifies stress in the chordwise direction, is the ratio of 
chordwise compression P x to the buckling value of chordwise stress 
which corresponds to zero stress in the spanwise direction. Positive 
values of this parameter indicate chordwise c empress ion and negative 
values imply tension. Curves are shown for two pinned-edge panels with 
aspect ratios of 0.5 and 1. 
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It must be remembered that the calculations apply to flat panels 
only. Hence, the curves in figure 3 are valid only up to the chordvise 
load which results in buckling. If the midplane stress in the spanwise 
direction is zero, the curves are therefore valid, by definition, only 
up to the compression parameter equal to one. If tension is applied to 
the panel in the spanwise direction, the chordwise load necessary to 
cause buckling is raised j the curves would then be valid to a higher 
value of the compression parameter. Spanwise compression would have 
the opposite effect. It should be noted that this is the only influence 
of spanwise stress on the flutter boundaries of flat panels. (See 
ref . 26 . ) 

The influence of chordwise stress is shown by the variation of the 
curves themselves. Note especially that, for each aspect ratio, there 
is a critical value of chordwise compression for which the theory requires 
very large thicknesses to prevent flutter. In real panels, the required 
thickness in this range is dependent on the amount of structural damping. 
However, it is apparent that these critical combinations of aspect ratio 
and chordwise compression should be avoided in design. 

The mode shapes of flutter of flat panels have an interesting feature 
that is illustrated in figure 4. Calculated chordwise variations of typ- 
ical flutter-mode shapes are shown for panels with both pinned and clamped 
edges. The air flow is from left to right. Note that the motion is con- 
centrated toward the rear of the panel for both panels. This "tail- 
wagging" type of motion is characteristic of the observed motion in actual 
instances of panel flutter. 


BUCKLED PANE IS 


Consideration will now be given to the flutter of panels which have 
been buckled by heating or by the application of midplane loads . Buckled 
panels have been treated in a number of theoretical investigations, most 
of which are confined to panels of infinite aspect ratio. (See refs. 1, 

2, 3, 9, 11, 23 , and 26 . ) In addition, a few experimental studies have 
been made. (See refs. 13, l 6 , and 23 .) 

The effect of buckling may be seen from the results in figure 5 - 
The required thickness ratio t/a is plotted against Mach number for 
steel panels of infinite span in air at sea level. Curves are shown at 
high Mach numbers for both pinned and clamped panels. The upper pair of 
curves (fig. 5 ) are calculated boundaries which apply to buckled panels 
of infinite span without regard to buckle depth. They result from closed- 
form solutions based on the so-called "transtability" concept. (See 
ref. 1.) The lower pair of curves (fig. 5 ) are the corresponding results 
for flat panels. Note that the buckled panels require much larger thick- 
nesses to prevent flutter than do the unbuckled panels. 
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Experimental results for clamped, buckled panels are also shown in 
figure 5 . With the exception of the point at M = 2.18, these results 
were taken from reference 16. (They have been corrected to apply to 
steel at sea level.) The point at M = 2.18, which was obtained inde- 
pendently (ref. 23 ), appears to confirm the other experimental results. 
The dashed theoretical curve for clamped buckled panels differs on the 
conservative side by about 12 percent. 

In references 13 and 16, some experimental results were presented 
for a group of clamped rectangular panels buckled by heating. These 
results are repeated in figure 6 . 

Figure 6 is constructed with modified -thickness -ratio parameters 
for both the ordinate and abscissa. One thickness ratio is based on the 
chord and the other, on the span. The data are applicable to all Mach 
numbers greater than about 1.3* The construction of the plot is such 
that data for panels of a given aspect ratio fall on a radial line from 
the origin. The aspect ratios and orientation of the test panels are 
shown by the small figures at the outer ends of the radial lines. In 
each case, node lines are shown to illustrate the observed buckle pat- 
terns. The distance outward from the origin of the plot is a measure 
of the panel thickness. Solid points designate panels that fluttered 
and open points indicate that no flutter was observed. 

An attempt has been made to calculate a theoretical flutter boundary 
for this particular group of panels buckled by heating. Because the 
midplane stresses in the test panels had not been measured, it was neces- 
sary to estimate these stresses from the observed buckle patterns. Some 
details of this calculation are given in the appendix. The solid curve 
is the resulting approximate boundary. The curve has been shown dashed 
on the lower right because the calculation is a two-mode approximation 
which must be considered unreliable for the very low-aspect-ratio panels. 

The agreement in figure 6 between the theoretical boundary and the 
experimental data encourages the viewpoint that successful analyses of 
rectangular buckled panels are possible. Of particular interest is the 
finding that, as in the case of flat rectangular panels, there are crit- 
ical combinations of aspect ratio and stress for which even thick panels 
can flutter. One such critical point occurs in the calculated boundary 
near an aspect ratio of 0.5* Perhaps this explains why the test panels 
with an aspect ratio of 0.5 were more prone to flutter than the other 
test panels . 

It must be pointed out that the theoretical boundary is quite sensi- 
tive to the type of buckling load, especially for aspect ratios between 
0.5 and 2. Consequently, this boundary applies only to the particular 
series of tests shown. It should be emphasized that this boundary is not 
intended to be a universal flutter boundary for finite-aspect-ratio 
buckled panels. In fact, no such universal boundary exists. 
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CIRCULAR CYLINDERS 


One other configuration which has received some preliminary atten- 
tion is the infinitely long thin-walled circular cylinder with axial 
flow over its outer surface. (See refs. 14, 21, and 27.) Typical 
results , for empty steel cylinders, are shown in figure J. 

The ordinate in figure 7 is the ratio of wall thickness to cylinder 
radius t/r and the abscissa is the Mach number. The curves are theo- 
retical flutter boundaries which apply to unstiffened infinitely long 
cylinders. The solid curves are boundaries that result when there is 
no static-pressure differential across the wall; the dashed curves, on 
the other hand, correspond to the existence of enough internal pressure 
to cause a circumferential tension near the ultimate. 

It will be noted that the required wall thicknesses of unstiffened 
infinitely long cylinders increase rapidly with increase in Mach number 
and very quickly reach prohibitive values, even at very high altitudes. 
Internal pressure is seen to produce only small percentage reductions 
in these large thicknesses. Furthermore, it has been pointed out in 
reference 27 that negligible benefit can be expected if internal damping 
is taken into account. 

On the other hand, an effect that shows promise of reducing the 
thicknesses to reasonable values is the effect of finite length. For 
the purpose of illustrating this effect, an approximate calculation has 
been made for a steel cylinder with rigid ring-stiffeners spaced a radius 
apart. The calculation is based on equation (42) of reference 21 with 
the air forces replaced, for simplicity, by their plane static approxi- 
mation. This simplification imposes the requirements that the Mach 

number M be large and that M — » ^- where n is the number of full 
° a 2 it 

waves around the circumference of the cylinder (corresponding to the 
maximum required thickness for the prevention of flutter) and a is the 
distance between stiffeners. The resulting thickness ratio required to 
prevent flutter at M = 6 of the stiffened steel cylinder at an altitude 
of 35*000 feet with no internal pressure is shown by the small circle in 
figure 7 . This result is only approximate because the above conditions 
are only approximately satisfied (maximum t/R corresponds to n = 26 ) . 
However, it has the correct order of magnitude. Furthermore, in contrast 
to the extremely large thickness ratios predicted at M = 6 for the 
unstiffened infinitely long cylinder, the result for the stiffened cyl- 
inder has a reasonable order of magnitude. Thus, the large thicknesses 
predicted by analyses of infinitely long unstiffened cylinders apparently 
do not apply to practical configurations of finite length. 


CONFIDENTIAL 



8 


CONFIDENTIAL 


NACA PM L57D24c 


CONCLUDING REMARKS 


It may be stated that much progress has been made in the solution 
and understanding of the panel flutter problem. Although, in general, 
panel flutter is of concern only in the design of very thin panels, there 
appear to be critical combinations of panel aspect ratio and midplane 
stress for which even very thick panels may flutter. It should be 
emphasized that there is a need for more work on certain phases of the 
panel flutter problem. In particular, there is a need for further con- 
sideration of cylinders of finite length and for more experimental 
results for assessing the validity of theoretical analyses. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va . , March 7 > 1957 • 
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APPENDIX 


TRANSTABILITY ANALYSIS OF THE FLUTTER OF BUCKLED 
RECTANGULAR PANELS WITH CLAMPED EDGES 


The application of the "transtability" concept (ref. l) to rectan- 
gular buckled panels with pinned edges is discussed in reference 26 . 

The corresponding results, for clamped panels, are presented in this 
appendix . 

' The transtability problem is governed by the plate equation 


DV W + P x w xx + P y w yy 


2q 
Jm 2 - 


w x = ° 


( 1 ) 


where w is the panel deflection and the subscripts denote differentia- 
tion, and by the proper boundary conditions. The assumption is made that 
the panel deflection shape is adequately represented in the chordwise 
direction by a linear combination of the first two modes of a uniform 
clamped-clamped beam and in the spanwise direction by a single mode, the 
pth mode. Then, by the Galerkin method, the condition for the existence 
of a nontrivial solution is found to be 

iM* 2 - (4i + Vn + b p) (4' + v*? + ®p) (2) 


where 


A = 2qa 


Jm 2 - 1 D 


P x a 2 

Ap d +2 ' 

/afi (2) 
VhJ T>P 

w (A 2 T {2) 

Bp ~ (t j X PP 

P y &2 / a \ 2 

D + (5) % 
_ - L pp . 

4l> - £ 
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The function 0 m (i) is the mth mode for a beam of unit length. The 

beam modes and the integrals are given conveniently in refer- 

ences 29 and 30 . 


As pointed out in reference 2 6 , the critical value of the dynamic- 
pressure parameter A depends on the midplane loads P x and Py applied 

to the panel at the boundaries or induced by heating a panel with fixed 
boundaries. Let Py = RP X . Substituting Ap and Bp from equations ( 3 ) 

into equation (2) yields a quadratic equation for the buckling loads P x 

corresponding to given values of R and A. As A increases from zero, 
the buckling loads approach each other until at the critical "transtability" 
value of A they coalesce and disappear. This condition is given by the 
vanishing of the discriminant of the quadratic which yields the result 



Calculations have been made of the quantities: 


and 



for a series of test panels buckled by heating. The integer p, speci- 
fying the number of half-waves in the spanwise direction, was taken as 
the number of spanwise half-waves in the observed buckle patterns. 
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Furthermore, since the actual loads in the test panels had not been 

P 

measured, the ratio R = — — was estimated from the buckle patterns 

P 

with the aid of figure 5 of reference 31. Specifically, for each aspect 
ratio, the average longitudinal buckling load corresponding to the number 
of observed spanwise or chordwise half-waves was used along with the 
corresponding transverse load. For an aspect ratio of 1, it was assumed 
that R = 1 (P x = P y ) . The resulting calculated boundary is shown in 

figure 6 along with the test results. 
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REQUIRED THICKNESS OF INFINITE -ASPECT-RATIO 
UNSTRESSED STEEL PANELS AT SEA LEVEL 



EFFECT OF ASPECT RATIO ON REQUIRED THICKNESS OF 
UNSTRESSED SIMPLY-SUPPORTED STEEL PANELS 
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EFFECT OF MIDPLANE STRESS ON FLUTTER 
OF UNBUCKLED PANELS 




(BUCKLING Px) Py=0 


Figure 3 


TYPICAL FLUTTER -MODE SHAPES 



LEADING EDGE TRAILING 


PINNED 


CLAMPED 


EDGE 


Figure 4 
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EFFECT OF BUCKLING ON REQUIRED THICKNESS OF 
INFINITE- ASPECT-RATIO STEEL PANELS AT SEA LEVEL 



REQUIRED THICKNESSES OF A SERIES OF 
TEST PANELS BUCKLED BY HEATING 



Figure 6 
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REQUIRED THICKNESS OF INFINITELY 
LONG STEEL CYLINDERS 



Figure 7 
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